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Abstract 



We discuss the Gupta-Bleuler quantization of the free electromagnetic 
field outside static black holes in the Boulware vacuum. We use a gauge which 
reduces to the Feynman gauge in Minkowski spacetime. We also discuss its 
relation with gauges used previously. Then we apply the low-energy sector of 
this field theory to investigate some low-energy phenomena. First, we discuss 
the response rate of a static charge outside the Schwarzschild black hole in four 
dimensions. Next, motivated by string physics, we compute the absorption 
cross sections of low-energy plane waves for the Schwarzschild and extreme 
Reissner-Nordstrom black holes in arbitrary dimensions higher than three. 
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I. INTRODUCTION 



The prediction that black holes should thermally evaporate has sparked much interest 
in quantum field theory in curved spacetimes. One of the difficulties in studying fields in 
Schwarzschild and other black hole spacetimes, even when the fields are non-interacting, 
stems from the fact that the solutions to the field equations are functions whose properties 
are not well known.Q In the low-frequency regime, however, the situation is much simpler in 
Schwarzschild spacetime. In this regime, the mode functions of the massless scalar field are 
well known Q . The present authors and Sudarsky used this fact to find the response rate of 
a static scalar source |^ and an analytic approximation for the emission rate of low-energy 
particles from classical sources outside the Schwarzschild black hole in closed form. On the 
other hand, the field equations of the electromagnetic field in a black hole spacetime are not 
decoupled and are difficult to analyze in the Lorenz gauge. However, if we require the field 
to be divergence- free on a two-sphere (the spherical Coulomb gauge), the equations for the 
physical modes reduce to decoupled scalar field equations. Furthermore, solutions in terms of 
familiar special functions can be found in the low-energy regime. These observations enabled 
us recently to calculate the response rate of a static electric charge outside a Schwarzschild 
black hole in closed form 0. 

In this paper we examine free quantum electrodynamics in static spherically symmetric 
spacetimes of arbitrary dimensions in a modified Feynman gauge. (This gauge is closely 
related to the Aq = gauge used by Cognola and Lecca ||^ and reduces to the Feynman gauge 
in Minkowski spacetime.) Then, we calculate some low energy quantities in electrodynamics 
outside spherically symmetric black holes. First, we review the calculation of the response 
rate of a static charge outside the four- dimensional Schwarzschild black hole in the Unruh 
vacuum 0. Next we calculate the low energy absorption cross sections of photons for the 
Schwarzschild and extreme Reissner-Nordstrom black holes in arbitrary dimensions higher 



than three, extending some results obtained by Gubser ||T0[ using a method |[TT|Jl^ based on 



the Newman-Penrose formalism [|I3 



The paper is organized as follows. In Section |I| we present the mode functions of the 
electromagnetic field in the spacetime of a spherically symmetric black hole in our modified 
Feynman gauge. Then we discuss the corresponding quantum theory and show how the 
Gupta-Bleuler condition (see, e.g., |T^) is implemented to obtain the physical states. In 
Section |ITI| we compare the physical modes in the spherical Coulomb gauge [0 with the 



ones obtained in the modified Feynman gauge. In Section |I^ we review the calculation of 
the response rate of a static charge outside a four- dimensional Schwarzschild black hole. In 
Sections and we present the photon absorption cross sections by the Schwarzschild 



and extreme Reissner-Nordstrom black holes of arbitrary dimensions higher than three. In 
Sec. [VII| we summarize the main results and make some remarks. In Appendix ^ we compute 



some components of field-strength two-point function in Minkowski spacetime using spherical 
polar coordinates and show that they agree with those obtained using Cartesian coordinates. 
In Appendix ^ a summation formula for Legendre functions used in Section |V] is derived. 
In Appendix |C|, a formula which relates the absorption probability to the absorption cross 



^See Ref. for some known properties in the Schwarzschild case. 



2 



section is derived in arbitrary dimensions. We use the metric signature (H . . . and 

the natural units with G = h = c = 1 throughout this paper. 



II. GUPTA-BLEULER QUANTIZATION IN A MODIFIED FEYNMAN GAUGE 



In this section we analyze the field equations for the electromagnetic field in spherically 
symmetric and static spacetimes in a modified Feynman gauge. Then we discuss the Gupta- 
Bleuler quantization in this gauge. 

The line element of the spacetime we study is 



c/r^ = f{r)dt^ - h{r)dr'^ - r^dsl 



(2.1) 



where dSp is the line element of a unit p-sphere. We assume that /(r) and h{r) are positive 
for r > th and that both /(r) and h{r)~^ have simple zeroes or both have double zeroes 
at r = Th, where th is the horizon radius. We also assume that f{r),h{r)^^ — > 1 as 
r — s> oo. (Most results in this section, however, are independent of these assumptions.) Let 
us introduce the Wheeler tortoise coordinate r* by 



dr* 
dr 



(2.2) 



Then r*(r) is a monotonic function with domain {th, +oo) and range in (— oo, +cxd). Let us 
define for any two functions qi{r*) and q2{r*) the inner product 



ill, 12) 



dr*qi{r*)q2{r*) 



(2.3) 



where the overline denotes complex conjugation. 

The Lagrangian density of the electromagnetic field in a modified Feynman gauge is 



-9 



--F.^F^" - -G^ 
4 ^ 2 . 



with 



G = VM^ + K^A^ , 
where the vector is independent of A^. Hence the equations of motion are 

V^F*^^' + V^G - K^G = . 

Here we choose 

^^ = (o,/7W,o,o), 

in which case Eq. ( p.5|) is written as 



G 



dr 

hrP 



fh 



Ar 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 
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Here i denotes angular variables on the unit p-sphere with metric fjij and inverse metric 
fj^^ [with signature (+ ■ ■ ■ +)], Vj is the associated covariant derivative on 5*^ and V* = fj^^Vj. 
This choice for K'^ is convenient because the equation for At decouples from the other ones. 
The field equations (|2.6| ) become 




+ 4v' A = , 



1 - 



f 



/ 



— dfAi+ ^ 



/ 



(2.9) 

(2.10) 



(2.11) 



where = f/ijV^V-'. 

We shall now describe a complete set of solutions We assign A the value 

for what we call the non-physical modes, 1 or 2 for the physical modes and 3 for the pure- 
gauge modes. (These modes will be given below.) The label n distinguishes between modes 
incoming from the past null infinity (denoted with n =^) and those coming out from 
the past horizon (denoted with n =^).R The solutions with ^(0"''^''^) = 0, (/i ^ t), and 



A 



(On; uilm) 
t 



-iuit 



(2.12) 



will be called "non-physical modes" because they satisfy the field equations ( p.9|) -( pTTD but 
not the gauge condition G = 0. Here Yi^ is a scalar spherical harmonic on the unit p- 
sphere with V^y^m = —li} + P — 1)^;™, where / = 0, 1, 2, . . . and m denotes a set of p — 1 
integers (mi, . . . , mp_i) satisfying / > mp_i > . . . > m2 > \mi\. (See Ref. for a concise 
description of spherical harmonics on the p- sphere.) They are normalized as 



/ 



dfipYimYi'. 



(2.13) 



where d^p is the volume element of the unit p-sphere. The function R^^P\r) satisfies 
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+ 



h rP dr 



d 



l{l+p 



fh dr 



p(On) 



(2.14) 



We will determine the normalization of the functions R^i^\r) later. The pure-gauge modes 
are given as 



(2.15) 



^Here we treat only the solutions proportional to e with u> ^ 0. Thus, if there is any nonzero 
static field, we will be considering fluctuation about that solution. 
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where 



(2.16) 



(As usual, the pure-gauge modes satisfy the field equations (pl9| ) - (|2.11|) and G = 0.) The 
other independent solutions (A = 1,2), which represent physical degrees of freedom, i.e., 
which satisfy the field equations and G = but are not pure gauge, will be chosen to have 
At = 0. Those with A,, are given as 



where 



'uj^ _ l{l+p-l] 
J ~' 



(In) 



V) + — — 

^2 



h dr \\/fh 



^ ~(ln). 



0. 



(2.17) 



(2.18) 



Note here that the condition G = cannot be solved if / = 0. Hence, we have Z > 1. The 
corresponding angular components can be found by solving the condition G = as 



A 



{In; ujlm) 



l{l + p-l)\l hdr 



(2.19) 



We call these modes "physical modes I" . The other set of physical solutions can be obtained 
by letting At = A,. = and 



A 



(2n; ujlm) 



p(2n). . 



where 



UJ 

1 



+ 



d 



' fh rP~2 dr \\ h 



dr 



i?2")(r) = 0. 



(2.20) 



(2.21) 



Here, the Y^''"^^ are divergence-free vector spherical harmonics on the unit p-sphere satisfying 



V^VkYl'""^ - V^ri'"^) = -(/ + l){l+p- 2)yI'"'^ (2.22) 



and 



(2.23) 



(See, e.g., Refs. [|T^,|T^].) We call these modes "physical modes H". Physical modes I and 
II obtained here (and restricted to four dimensions) are identical with those in the Aq = 
gauge g]. 

In order to discuss Gupta-Bleuler quantization of this field it is convenient to introduce 
a generalized Klein-Gordon product of classical solutions of Eq. (|2.6|) . We first define 



H 



^d[V,A, 



(2.24) 
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where G is given in Eq. (|2.5|) . Note that \/—g H*'^ is the canonical conjugate momentum of 
Ay. We write 

n(c)/^. ^ n^1^^=4c) (2.25) 

for any solution A^^^ = ^{An;wim)_ p^^, (complex) solutions A'^p and A^^'-* we define 

W^[A^^\ = z [SPn^^')^'^ - im^Ap] . (2.26) 



The field equations ensure that this current is conserved. As a result, the generalized Klein- 
Gordon inner product defined by 



(2.27) 



where dS^ = daUfj,, is independent of the Cauchy surface S |T^. (Here da is the volume 
element of the Cauchy surface S with a normal unit vector n^.) Note that (A^'^\ A^'^ ^) = 
(y4(^'), y4(^)). This guarantees that the norm defined through ( |2.27D is real (and positive 
definite for the subset of physical solutions with positive frequency). By working explicitly 
with the definition (|2.27| ), we find in general that on a t = const surface 



(2.28) 



It is important to note that pure-gauge modes are orthogonal to any mode satisfying 
G = and, as a result, V^F'^^ = 0. This can be shown as follows. Suppose that A^ and 



^(f') satisfy the condition G = 0. Then since H^^)'^^ = - F>"'\ 



and similarly for 



the inner product ( p.27| ) can be written as (A^), A(f')) = (A«), A«'))inv where 



Now, let A^ = V^A'^'') be a pure-gauge mode. Then 



(2.29) 



(2.30) 



since F^'^">'^^ is ant i- symmetric. Q Thus, {A'^'^\ A'^'^'^) = if A^^-* is a pure-gauge mode and 
is a physical or pure-gauge mode, i.e.. 



^yj^{3n; wim) ^(3n'; Lj'/'m') ^ g 



(2.31) 



^j^(yin\ u)lm) ^{In'; oj'/'m') ^ g ^j^{3n; Lolm) ^(2n'; oj'i'm') ^ g 



(2.32) 



^Note that Eq. (|23o| ) would hold even if A(0 was an arbitrary function. This shows that 
inv is gS'Uge invariant. 
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Next let us examine the inner product of the non-physical solutions. By letting 

i?ir(r)-A/7r-/Vir(r), 



we find from Eq. ( |2.14| ) 



(On) 



dr* 



where r* is given in Eq. ( p.2|) and 



Vo[r*ir)] = f- 



h 



- Voir*) 



(2.33) 



(2.34) 



4^2 



2/2 4r V / /i / 2/ Afh 



For spacetimes where /(r) and h{r)~^ have simple zeroes at r = th, we find Vo[r*(r/i')] ^ 
in general. If Vo[r*(rj^)] > 0, then the non-physical modes coming out from the horizon have 
frequencies satisfying u"^ > Vo[r*(r//)]. In particular, for Reissner- Nordstrom spacetime in 
(p+2)-dimensions [|l^], where 

fir) = h{r)-' 



1_ 

r , 



r 



with r^^ = M ± y/M'^ - Q'^ (so that vh 



Volr'-irn)] 



(2r+)^ 



we find 



1- - 



p-1 



As expected, we have {f^^i^\f^3T^) = O5 i-^- the solutions of Eq. ( p.34| ) incoming from 
H~ and those incoming from J'~ are orthogonal to one another with respect to the inner 
product defined by Eq. (|2.3| ). By normalizing these solutions so that 



(0-*) 



I ~ ( iCur* I T) (0^) —iuir* 

uj/u) e +nii 'e 



(r' 
+00) 



-00) 



(2.35) 
(2.36) 



up to a phase factor, where uo'^ = uo'^ — Vo[r*(rj:/)] > and TC^i and TC^i are constants 
(with |7^S'^)| = |7^i^)|),wehave 



(^r\^ir) = (v^r\v^ir)=2^5(---o. 

Now, by using Eq. ( p. 2^ ) we find 

= -47ru;5„„/(5ii'5mm'5(t^ - to') . 



(2.37) 



(2.38) 



where we have used 5{uj - u') = {lj/uj)5{uj - cu') . Noting that A^^";'^'"') and ^(^"^^M have 
the same t-components, we immediately obtain from Eqs. ( p.28|) and (|2.31|) that 
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^j^iOn-Mm)^ ^(3n';.'Z'm')) ^ _4vr^5„„,5;;,<5^„,<5(cU - Uo') . (2.39) 

Next let us examine the physical modes. By letting 



i?i-)(r) ^ ^I^IZ^ (//.)^/V-/-Vir (r) , (2.40) 

UJ 



we find from Eq. (^]T 



(-^ + |^-V^1(0)^S"^ = 0, (2.41) 

where 

, J(/+p-l) p(p-2)f (p-2) f ff h'\ 

(We note that Vi — * as r ^ r+ for Reissner-Nordstrom spacetime unlike Vq.) We first note 
that the modes ^(i";^'™-) are orthogonal to both the non-physical and pure-gauge modes. 
Next we note 

(^(in;.M, Ai'nwi'm'y^ = • (2.43) 
Thus, by normalizing <^^^"''(r) as 

^ap) ^ ^^.r' ^ 7^(;p)e— * (r* ^ -oo) , (2.44) 

^a-) ^ g-...' + 7^i^)e-'■* (r* ^ +00) , (2.45) 

we have 

^j^{ln-Mm)^ ^ ^T,^Snn'5lv5mm'5{u - u') . (2.46) 

Finally, we let 

i?if(r)-r-(-^)/Vif^(r). 



Then from Eq. (|2.21|) , we have 



'^' + :£^-V^2(r*)V2") = 0, (2.47) 



dr 



where 



(We again note that V2 as r — ^ r+ in Reissner-Nordstrom spacetime.) Notice that 
Vi = V2 li p = 2. Hence, in the four-dimensional case (pl^P\r) and (p^^P\r) satisfy the 
same equation. These modes can easily be shown to be orthogonal to the modes previously 
discussed. By normalizing the functions f^^P^ {r) in the same way as f^^P'' (r) [see Eqs. ( p.44|) - 
( p.45| )], we have 
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(2.49) 



In order to quantize the field A^, we impose the equal-time commutation relations on 
the field and momentum 11*^ operators: 



i^(t,x),i,(t,x')l = ffl*'^(t,x),fl*'^(t,x' 



0. 



(2.50) 



A^(t,x),ff'^(t,x'^ 



-^5^'+Vx-x'l 



-9 



(2.51) 



where x and x' represent all spatial coordinates. The field can be expanded using the 
modes we have obtained before: 



^.(t,x) = E/_ 



°° \ A'n\uj\ 



(2.52) 



where A^^f^ is proportional to e"*"^*, A[^'^f^ = A 



(up) 



aj^p and p labels discrete 
quantum numbers. The commutation relations (|2.50|) - (|2.51| ) are equivalent to the follow- 
ing commutation relations in the "symplectically smeared" form as is the case for scalar 

fields m-. 



(A(^),i),(i,AK'))l =(AK),A(^')). 



(2.53) 



Since the inner product must be t independent, the inner product of A'-ZfP^ and A^^'p'^ can 
be nonzero only if u; = u'. Thus, we can write 



{A^^p\ A^^'p'^) = MPp'5{uj - uj') 
By using Eq. (p.54|) in (|2.53|) one finds 



^ 4t^ujMP'P'6{uj - u') 



Then, since Mpp is invertible in our case, we find 



O-ujp^ ^uj'p' 



(2.54) 



(2.55) 



(2.56) 



Note here that we immediately have [a^jp, a^/p/] = for oj.uj' > by letting oj' —uj' 
in Eq. ( |2.56D . By using the Klein-Gordon inner products computed above, we find the 
following commutators: 



{3n) (3n')t 



(In) {ln')t 



(On) (3n')t 
(2ra) (2n')t 



(2.57) 
(2.58) 



with all other commutators vanishing. The Gupta-Bleuler condition requires that any 
physical state |phys) satisfy 
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G^lphys) = 0, 



(2.59) 



where is the positive-frequency part of G = V^A^ + K^^A^j^. Since this quantity is 
nonvanishing only for ^(^On; t^im) ^ ^j^jg condition is equivalent to 

a2!;2|phys) = for all (n, u, I, m) (with u > 0) . (2.60) 

(We let u; > below.) The Boulware vacuum |0) ||20| is defined by requiring that it be 
annihilated by all a^^i^ operators (A = 0, 1, 2, 3). Note that the states obtained by applying 
any number of creation operators excluding a^^;^^ are all physical states. Any state of the 
form a^^;^^|phys) is unphysical because 



a, 



(On') (3n)t 



'u'l'm'(^llm Iphys) = -Snn'Sll'Smm'S{uJ - ^')|phys) ^ . 



Note also that the physical states of the form o^°m^lP^y^) have zero norm and are orthogonal 
to any physical states. Thus, as is well known, a physical state Iphys^) can be regarded as 
equivalent to any state of the form Iphys^) +a^°^^|phys2). We can take as the representative 
elements the states obtained by applying a^"^^^^, A = 1,2, on |0). 

Unphysical particles created by a^^i^^ will be in thermal equilibrium in the Hartle- 
Hawking vacuum ^T| for a static black hole if we require the gauge-fixed two-point function 
be non-singular on the horizons as in the scalar case p2| . There will also be a flux of unphysi- 



cal particles in the Unruh vacuum. Therefore, technically speaking, the Hartle-Hawking and 
Unruh vacua are unphysical if we impose the Gupta-Bleuler condition using the positive- 
frequency notion in the Boulware vacuum. Fortunately, the BRST quantization does not 
suffer from this problem since the physical state condition does not refer to any notion of 
positive frequency This issue, however, will not concern us in the following, since all 
our calculations will be at the tree level. 

Mode expansions in spherical polar coordinates are not very widely used and are quite 
different from the ordinary method using Cartesian coordinates. For this reason, we compute 
in Appendix A some components of field-strength two-point function in four- dimensional 
Minkowski spacetime using spherical polar coordinates. We find that they agree with the 
standard results obtained in Cartesian coordinates, as they should. 

III. THE PHYSICAL MODES IN THE SPHERICAL COULOMB GAUGE 

As we have seen, the modified Feynman gauge defined by Eq. ( pI7| ) is useful in finding 
explicit mode functions. We also noted that this gauge results in the same physical modes 
as in the Aq = gauge. There is another convenient gauge for finding physical modes, 
i.e., the spherical Coulomb gauge, V*v4j = 0. (One can readily show that any vector is 
gauge-equivalent to a vector satisfying this gauge condition.) Let us discuss the relation of 
the physical modes obtained in this gauge and those found in the previous section. R 



^ Again, we only treat solutions proportional to e with lo ^ 0. 
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It is clear that the physical modes ^(^2n;tjZm) ^^le previous section are also in the 
spherical Coulomb gauge. The other independent physical modes must have = in this 
gauge. First we note that if At and Ar are assumed to be spherically symmetric, then we 
find from the equations V^F^t = V^F^r = that Ftr is t-independent and {r^ / ^/Jh)Ftr 
is r-independent. Hence, if in addition the modes are assumed to be proportional to e"*"^* 
with a; 7^ 0, they must be pure gauge because Ftr = 0. (A time-independent solution has 
Ftr oc \/Jh/r^. This will represent a static Coulomb field, which we will disregard here.) To 
find the solutions which are not spherically symmetric we let At, Ar oc 11^6"*'^*. From the 
equation V^F^i = we find 



iujAt 



dr 



'IrP-'Ar 



(3.1) 



By substituting this and A = in V^F^,. = 0, we find that A^^'"'"^'™) = i?J,V"^(r)l^„e'^'^* 
(here we use primed indices to refer to the spherical Coulomb gauge), where 



_ l{l+p-l) 
f r"^ 



h rP~'^ dr \\ h ^ 



0. 



(3.2) 



The component A 



(I'n; u)lm) 



can be found from Eq. ( |3.1| ). (Then, one can readily show that 



V^F^t = 0.) Indeed one can verify that these modes are related to the ones in the modified 
Feynman gauge by a gauge transformation: 



(3.3) 



where 



^(Ira; u)lm) 



' f d ( 



l{l + p-l)^ hdr \y/fh 



R. 



(In) 



-iuit 



(3.4) 



This gives us 



(I'n; uilm) 



A 



2-p 



R. 



(In) 



^'^r'^ 1 ^(In) 



-ILOt 



-ibjt 



(3.5) 
(3.6) 



IV. RESPONSE RATE OF A STATIC CHARGE OUTSIDE A 
FOUR-DIMENSIONAL SCHWARZSCHILD BLACK HOLE 

Here we briefly review, in the context of the previous two sections, the calculation of 
response rate in the Unruh vacuum of a static electric charge in Schwarzschild spacetime 
performed in Ref. 

The line element of the four-dimensional Schwarzschild spacetime is 
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dr^ = f{r)de - f{r)-Mr^ - r^dsj (4.1) 

with /(r) = 1 — 2M/r, where ds^ = dO"^ + siv? dcjp' is the hne element of the unit two- 
sphere. (The horizon radius is th = 2M.) We compute the response rate of the static 
electric charge to the Hawking radiation present in the Unruh vacuum. The charge is placed 
at (r, 9, 0) = (ro, 6*0, 0o) as described by the current density 

= ^t-^A^ - - ^o)5(0 - 0o) . (4.2) 

sm U 

Since this current density is static, it couples to photons with zero energy. It turns out that 
the rate of response of this current to a single photon vanishes. However, the Bose-Einstein 
distribution of the thermal photons coming out of the horizon diverges at zero energy. This 
makes the rates of absorption and stimulated emission indefinite. This ambiguity can be 
resolved by replacing the current density ( ^4.2| ) by = (j*, j'", 0, 0), where 

f = "^r^'f' sir - ro)Sie - e,)5{<P - 0o) , (4.3) 
sm U 

f = .'^^/^ 9(r - ro)5i9 - 9o)Si^ - 0o) • (4.4) 

sm 9 

Here, 9(x) = 1 if x > and 9(a;) = if a; < 0. (The factor of is necessary to make 
the time average of the squared charge equal to g^.) We take the limit E ^ in the 
end, assuming that the rate is continuous at -E = 0. This continuity has been verified in 
Minkowski spacetime [Q. 

Since the current ( ^^ is conserved, it does not interact with pure-gauge particles corre- 



sponding to the creation operators a^^hn ■ It does interact with the states created by 
but these have zero norm and do not contribute to physical probabilities. Hence, as usual, 
we need to consider only the physical modes. Note also that only physical modes I are 
relevant because At = Ar = for physical modes II. 

Now, the transition amplitude from the Boulware vacuum |0) to the state \ln; ujlm) = 
a^^i2^\0) by the classical current density is given by 

^"zrn = ^ / d^x^filn; ujlm\A^\0) . (4.5) 

Let us recast Eq. ( [4.5| ) as 

A:^i^^2mr:,Jiuj-E), (4.6) 

where will be calculated later. The corresponding transition rate, i.e. transition prob- 
ability per asymptotic proper time, is expressed as 

T^Zin. = - E) ■ (4.7) 

Now we take into account that in the Unruh vacuum there is a thermal flux of photons 
coming out of the horizon with inverse temperature (3 = 8nM (with no particles coming 
from infinity). In the limit E ^ 0, only the absorption and stimulated emission of these 
thermal photons contribute to the response rate. For this reason, we treat only the modes 
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coming out of the horizon (n =— >)• The total response rate due to these modes with fixed 
angular momentum, in which absorption, spontaneous and stimulated emissions are all taken 



into account, is 



Lolm 



RElm - 

Thus, in the limit ^ we find 



+ 1 ) = 27rcoth^|rr 



Elm\ 



Elm\ 



V(2ME) . 



(4.8) 



(4.9) 



Next will be calculated. The equation satisfied by the function q^i{r) = 

{2M/r)(f^^i^\r) can be obtained from Eq. (|2.41| ) with p = 2 as 



d_ 

dz 



U}1 



dz 



+ 



1(1 + 1] 



z + 1 



z-1 



C = 0: 



(4.10) 



where z = r/M — 1. For small u, the wave is almost completely refiected back to the horizon. 
This implies that 



-2 sin[co'(r* — r^)] ^ —2Ljr* + const. {r — Th <^ ^^'"h; \ujr*\ <^ 1) , (4.11) 



where are constants. The minus sign has been inserted for later convenience. The Wheeler 
tortoise coordinate r* defined by Eq. ( ^^ is written in Schwarzschild spacetime as 



r + 2M In 



/ r 



\2M 



1 



Thus, Eq. ( 4.11 ) can be written as 

-4:Muj\n{z - 1) 



r — Th 't:uj rj^, \ujr*\ <^ 1) . 



(4.12) 



(4.13) 



Now, Eq. (^4.101) can be solved explicitly for u = 0. By choosing the solution for q^i such 
that v^^^P"* = [r / {2M)]q^i tends to zero as 2; — +00, we find in the small uu limit 



4:Mu{z + 1) 



iz~l)dQi{z) 
+ dz 



(4.14) 



where this has been normalized such that Eqs. ( |4.13 ) and ( ^4.14 ) are in agreement at r ^ rj:^. 

From gauge invariance of the amplitude ( [4.5|) it is clear that we may use the modes in 
the spherical Coulomb gauge, A^^ ~*\i^im)^ place of those in the modified Feynman gauge 
/[O-^'^^'^)^ Then the contribution from the r-components will be suppressed in the low 
energy limit due to extra factors of u [see Eq. (|3^)]. Therefore we need to consider only 
the t-component in this limit. The t-component can be found in the small u limit from Eqs. 
(g), (^ and (^J§ as 



A 



(l'^;a;im) 



4:iuj{z - 1) dQi{z) 



dz 



Im \ 



-iwt 



(4.15) 
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where we have used the Legendre equation satisfied by Qi {z) . By using Eq. (|4.15| ) in Eq. 
and comparing the result with Eq. (|4.6| ), we obtain 



r 



2ujq{zQ - 1) dQi{zo) 



Lulm 



Y 



(4.16) 



^7il{l + 1) dzo 

where Zq = Tq/M — 1. Then we find the proper Tesponse rate of the charge as [see Eq. ([4.9p] 



Olm 



q\zo - 1) 



dzn 



\Y, 



(4.17) 



where the factor yf^ro) appears because here we are calculating the transition rate per 
proper time of the charge. We can sum over / and m by using 



E \YiU0o,<Po)\' = i2l + l)/{4n) 

m=—l 



(4.18) 



and 



^ 2/ + 1 
^ /(/ + 1) [ dz 



dQiiz) 



n 2 



2Qi(^) 

(^2 - 1)2 



(4.19) 



A derivation of this formula is given in Appendix B. Thus, the total transition probability 
per proper time of the charge is 



R 



tot 



0/m 



+ 00 I T3 

1=1 m=-l \l firo) 



27r2 



gi(ro/M-l) 



(4.20) 



where a(ro) = Mr'^'^/y firo) is the proper acceleration of the charge. 

It is of course possible to work with the modes in the modified Feynman gauge directly 
and obtain the same result. If we do so, the contribution will come from the r-component 
of the current because the t-component of physical modes I vanish. 



V. LOW ENERGY ABSORPTION CROSS SECTION OF THE 
SCHWARZSCHILD BLACK HOLE 

Here we consider the Schwarzschild black hole in j9 + 2 dimensions given by the line 
element ( f^.l| ) with /(r) = 1 — (r///r)^~^. If the absorption probability of the black hole 
associated with physical modes I and II are Vi^^ and "P^^ respectively, then the absorption 
cross section is given by 

-=;^f (5.1) 

where / and uj are the angular momentum and frequency, respectively^ of the modes and 
Vlp = 27T^'P^^^^'^ /r[{p + l)/2] is the volume of the unit p-sphere. Here, M; and m/^"* are the 
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the p-sphere. These multiphcities are given by (see, e.g., Refs. [15,16 



multiphcities of the scalar and divergence-free vector spherical harmonics, respectively, on 

(5.2) 



M, 



(1) 



{2l+p-l){l + p-2) 



and 



M, 



(2) 



{p-l)W. 

{2l+p-l){l+p-l)\ 
(/ + l)(/+p-2)(p-2)!(/-l)!- 



(5.3) 



Eq. ( p. ID is derived in Appendix C. 

For either physical modes I or II, we have for large r*) that [see Eqs. ( |2.41|) and 



l{l+p-l) + 



p{p - 2) 



(5.4) 



The asymptotic solution with the correct normalization, i.e. with the coefficient of e 
being one, is 



^ur\'r) ~ V27ro;rJ,+(p_i)/2(a;r) (r// < r) 



Thus, for small ur 



where 



Ci 



4:71 furny^P/'^ 



r(/ + (p+l)/2) V 2 



) 



(5.5) 



(5.6) 



(5.7) 



Now, near the horizon, where Vx <ti uj^ (A = 1,2), these solutions behave like [see 
Eqs. ([2:411 ) and (H^)] 



(5.8) 



where the constants Df''' will be determined later. Hence, for |ct;r*| ^ 1 we have in this 



region ipi^i '(r) d\'^^ and the absorption probability is given by Vl"^' = \Dl 

In order to find the coefficient dI^\ we solve the equations for v^^^^'*(r) with u = 0. For 
physical modes I, it turns out that the function r^R^^^'' satisfies a hypergeometric equation 
with the variable w = 1 — {r/ruY'^- Thus, we find in the small uj limit 

V^r\r) = d\'^ [rlmf-Pl'' F{~{1 +p- l)/ip - 1), l/ip -l);l;w). (5.9) 

This function approaches D;^'* for r — > r^/ as required. By using the formula |2^ 

r(7)r(/3-a) 



,(A) 



F{a,p;-f;w) 



r(/?)r(7 - a 



i-wy {a</3, -u;>l) 



(5.10) 
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in Eq. ( |5.9| ), we find asymptotically 



r/ru) 



l+p/2 



[th <^r <^uj 



where 



(1) 



By comparing Eqs. 



r[(2/ + p-i)/(p-i)] 
r[//(p-i)]r[(/ + 2p-2)/(p-i)] 

and ( ^.111) , we obtain 



V, 



(1) 



(1)12 



(5.11) 



(5.12) 



(5.13) 



Thus, the absorption probability behaves like 0;^'+^. 

Similarly, for physical modes II, we find in the small uj limit 

^T\r) = (r/r^)(^/^-^) F(-(/ + l)/{p - 1), (/ + p - 2)/(p - 1); 1; ^) . 

Then by using Eq. ( |5.10| ) we find asymptotically 



l+p/2 



where 



E-, 



(2) 



r[(2/ + p-i)/(p-i)] 



V[{l+p-2)/{p-l)]V[{l+p)l{p-l)]- 
By comparing Eqs. (5^) and ( 5.15|) , we obtain 



= [Q/E\ 



{2)i2 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



Again, the absorption probability behaves like lo^^^'^ . Therefore the physical modes with 
/ = 1 give the dominant contribution to the absorption probability, as expected. 

By substituting v'^^' and v'f^ in Eq. (|0| ), we obtain the total absorption cross section 
for low energy as 



a 



2{p+l) 



2p[T{p/{p-l))]' 
[T{{p+l)/{p-l))y 



+ 1 



(5.18) 



where Ah = ^pT^ is the area of the black hole. 

One would need to solve the differential equations satisfied by v'^^^'' numerically to find 
the low-energy absorption cross section for non-extreme Reissner-Nordstrom black holes.0 
However, it is possible to calculate it analytically for the extreme Reissner-Nordstrom black 
hole. We turn to this problem in the next section. 



^It is possible to obtain a result in a closed form if one takes the limit u ^ and r-| 



simultaneously |10| 
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VI. LOW ENERGY ABSORPTION CROSS SECTION OF THE EXTREME 
REISSNER-NORDSTROM BLACK HOLE 



The line element of the extreme Reissner-Nordstrom spacetime is given by Eq. 
with /(r) = h{r)^^ = [1 — (rj^/r)^^-*^]^. As in the Schwarzschild case, we have f^i^i*~\r) ~ 
Ci{r /rH)^^^^"^ for rn ^ r <^ u^^ where Ci is given by Eq. (|5.7| ). 

Now, let us analyze the physical modes close to the horizon. Concerning mode I, we 
write 



(In) 



(-a;V*^ > th 



(6.1) 



where we have used that the Wheeler tortoise coordinate r* [see Eq. ( |2.2| )] behaves in this 
region as 



Hence we obtain 



jp- 1) rj^ 
r ~ Th 



(6.2) 



D 



(6.3) 



with 



1 Ijl+P-I) 

4 {p-iy 



1 / 

+ 



2 p-1 



(6.4) 



This solution behaves like -D/^'*e (up to a phase factor) very close to the horizon {—ujr* ^ 
1), where the effective potential Vi becomes negligible, as expected. Hence, the absorption 



probability is vj;^"^ = \Dl'''\'^. Note that, for small z with positive z/, one has 



)(1)|2 



(6.5) 



Hence, 



{p - l)(r - rn) 



v-l/2 



,1/3 



< -cur* < 1 



(6.6) 



where 



(1) 



TT 



2(p-l) 



u-l/2 



(6.7) 



Again, in order to determine Dl , we solve the equation for ^plj^ {r) in the small uj limit: 



1 



(6.8) 
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where w = 1 — {r/rnY ^ i ^"^^ '^^ note that Eq. ( |6.8|) has been normahzed such that it agrees 
with Eq. ( |6.6|) near the horizon. On the other hand, we find asymptotically 



I 



l+p/2 



l+p-1 

By comparing this with Eq. ( |5.6D we find 

il+p-l)Q ' 



(6.9) 



-Pi 



(1) 



(1) 



(6.10) 



We note that Vi^^ oc cjP+^'p/^p Hence the dominant contribution to the scattering cross 
section comes from the / = 1 term: 



(1) 



2^2 



[{p - 1) V(P-i)r [{p + 3)/2] r [{p + l)/{2p - 2)f \ 2 J 



/cjrjy\2+p+2/(p-i) 



(6.11) 



The contribution from physical modes II can be found in a similar manner. By substituting 
Eq. ( |6.2D in Eq. ( ^^.47] ), we have near the horizon 



<f (/ + l)(/+p-2) 
dr 



*2 



[p — lyr 



2 /y* 2 



(2n) 



{^—uj'^r*^ ^ th 



Then ^p^^^'' will have the same form as >^^i 
replaced by 



in Eq. 



(6.12) 

except that v in Eq. ( |6.4| ) is 



I (/ + l)(/+p-2) 
4 (p-l)2 



nl/2 



For p > 3, the contribution from physical mode II is suppressed compared to that from 
physical modes I since v' > v. For p = 2, i.e. in four dimensions, Eqs. (|6.1| ) and (|6.12| ) 
satisfied by i^^^f' and ip'^f' are the same. Thus, the absorption probabilities of these two 
types of modes are equal for each /. 

Hence, we obtain the following results: 



(T = - {uthYAh , for p = 2 



(6.13) 



(p + 1) [2{p-l)] 



2/(p-l) 



[T[{p+l)/{2p-2)]] 



2 ' 



for p > 3 . 



(6.14) 



Eq. (|6.13| ) is in agreement with the result derived by Gubser |T0[ . It is interesting to note 
that the low energy absorption cross section is proportional to a fractional power of uj if 
p > 4. 
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VII. CONCLUSIONS 



In this paper we showed that the field equations for free electrodynamics can be re- 
duced to decoupled scalar field equations in a modified Feynman gauge in the spacetime 
of spherically symmetric black hole. (We also noted that the equations for the physical 
modes simplify in the spherical Coulomb gauge.) Then we examined the Gupta-Bleuler 
quantization in this modified Feynman gauge. 

Next we reviewed the calculation of the response rate of a static charge outside the four- 
dimensional Schwarzschild black hole. It is easy to extend our result to the Schwarzschild 
black hole in arbitrary dimensions, though we cannot simplify the resulting infinite series. A 
closed-form expression for the response rate of a static source has been obtained for massless 
scalar field in four- dimensional Reissner- Nordstrom black hole [^]. The extension of the 
results in this paper to the Reissner-Nordstrom black hole will also be possible in principle, 
but the result will probably not be expressible (even) as an infinite series of familiar special 
functions in any dimensions. 

Finally we calculated the absorption cross sections of low energy photons by the 
Schwarzschild and extreme Reissner-Nordstrom black holes in arbitrary dimensions. The 
corresponding cross section of massless scalar particles is known to be equal to the horizon 
area as long as the black hole is spherically symmetric. No such universality holds for 
photons as our results and previous results show. It is interesting that there are two 
modes, modes I and II, which behave differently in higher dimensions. It is also intrigu- 
ing that the absorption cross section scales as a fractional power of the energy a;, i.e. as 
^2+p+2/(p-i) p-|-2 dimensions. Finally, it would be interesting to compare the results here 



or similar results obtained using our gauge with those in string theory. (See, e.g., Ref. |10 
and references therein for examples of such comparison.) 
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APPENDIX A: FIELD-STRENGTH TWO-POINT FUNCTION IN SPHERICAL 
POLAR COORDINATES IN FOUR-DIMENSIONAL MINKOWSKI SPACETIME 



Minkowski spacetime is the simplest spherically symmetric spacetime. In this Appendix 
we compute some components of the two-point function 



G 



(Al) 



in spherical polar coordinates and verify that they agree with those previously computed 
in Cartesian coordinates in four-dimensional Minkowski spacetime. We need to consider 
only physical modes because pure-gauge modes do not contribute to F^j, and the coefficient 
operators of the non-physical modes have zero commutators with all operators appearing on 
the right-hand side of Eq. ( |A1| ). 

Since there is no horizon, all modes come in from infinity and go out to infinity. Therefore, 
we do not need here the label n =^ , — > which distinguished between the modes coming 
out of the horizon and those incoming from infinity. By the remark after Eq. ( |2.43| ) we 
have V'iV(^) — 2ti;rj;(u;r), where ji{x) is the spherical Bessel function of order I, in four- 
dimensional Minkowski spacetime. Hence for physical modes I we have exactly 



n , ^ 2.//(/ + l) 



(A2) 



where Yim{0, (j)) are the familiar scalar spherical harmonics on 5*^. The tr-component of the 
corresponding field strength is 



'tr ' = -I— J, 



ji{iur)Yirn{0,(j))e 



-iujt 



(A3) 



since the t-component vanishes. These modes are sufficient for computing Gtrtr{x, x') because 
physical modes II have Af = = 0. We obtain with x = (t, r, 9, cj)) and x' = {t', r', 9', 0') 



Gtrtr (x, X 



-| oo / 

1=1 m=—l 
oo 



irrr 



)Y, 



Im \ 



— ^Z(/ + l)(2Z + l)PKcos7)Q 



0') / dujujji(ujr)ji(iur')e-'''^'-''^ 
Jo 

-{t -t' - ief + + 



where we have used in the last step the formulas 



2rr' 



(A4) 



E 

m=—l 



)Y. 



Im \ 



2/ + 1 
47r 



P/(cos7) , 



with cos 7 = cos 9 cos 9' + sin 9 sin 9' cos (</ 
We will simplify Eq. (K3) later. 



and 25 



2rr 



-,Qi 



-it - f - ief + + r 
2rr' 
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Next we consider the component 



Ge,f,e<t>ix, x') 



sin 6 sin 6' 

4 



where the anti-symmetric tensor on 5^, e*-', has the component e^^ = (sin^^) ^. (This 
definition differs from that of Regge and Wheeler p8| by a minus sign.) Only the modes 
y^{2;ujim) contribute here. First, we have ^^^i{r) = 2ujrji{ujr). Then since R^^i{r) = ^l^i{r) 
for p = 2, we find 



^(2;.M ^ 2a;rj,(^r)r/'™V,0)e~'"*, 
where the transverse vector spherical harmonics y/'™^ are given by 



The corresponding field strength is given by 

l^,pi2Mm) ^ 2^l{l + l)urji{ur)YUeA)e 

Then in exactly the same way as for Gtrtr-, we find 
sin sin ^' °° 



(A5) 



(A6) 



-iujt 



Ge,j,e<i){x, x) 



5]/(/ + l)(2/ + l)P;(cos7)Qi 



'{t -t' - iey + + r 
2rr' 



(AT) 



(A8) 



On the other hand one can calculate Gap^s{x^ x') using plane wave modes. In the Carte- 
sian coordinate system this is given in components as (see, e.g., Ref. p9[| ) 



G^^Ux.x') = -47r{l + 2[(C)^ - {Ct?]r'}Gt{x,x') 



(A9) 



Gl2ix,x') = -8nUbr'Gtix,x') , 



(AlO) 



GlaL{^,x') = GiSi(x,x') = -87rCM-'G+{x,x') 



(All) 



G2b{^,x') = 4vr{l + 2[{Q' + {C,r]r'}Gt{x,x' 



(A12) 



G^f,l{x,x') = 87rUbr'Gtix,x') 



(A13) 



where a, b, c take values x, y, z with a, b and c being all distinct (no summation convention 
is used here) and = x°' — x". We have defined Gq{x, x') = (47r^^'^)^^ , where 



e = {t-t'- tef -{x- x'f - (y - y'f 
= [t — t — ie) — r — r + 2rr cos 7 . 



,z — z' 



i\2 



(A14) 
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The components Gtrtr in spherical polar coordinates can be obtained from the Cartesian 
ones (|A9| ) - (|A13 ) by using standard tensor transformation: 



Gtrtrix, x') = — - ({xx' + yy' + zz'\^ + 2 ixy' — x'yS^ + ixz' — x'zY 

+ [yz -yz) Jj , 

= I \{t -t' - ief - - r'^j cos 7 + 2rr'\ . 
7r"'4° L L J J 

Similarly, we find 



Gs60d,(x, x') 



2f6 



{ {t-t' - ief - - r'^ cos 7 + 2rr'} 



(A15) 



(A16) 



Eqs. ( A15 ) and and Eqs. ( A16 ) and (^) can be shown to agree by using the formula 

2{tz-l) 



Y^l{l + l){2l + l)Pi{t)Qi{z) 



(A17) 



which can be proved by applying {d/ dt)[{t'^ — l)d/ dt] on both sides of the well-known formula 



Y,{2l + l)Pi{t)Qi{z) = ll{z-t) 
1=0 

and using the Legendre equation satisfied by Pi(t). 



(A18) 



APPENDIX B: A DERIVATION OF EQ. ( gTTQl) 



Define 



°° 2/ -I- 1 



By the Legendre equation satisfied by Qi{z), we have 

1 2' 



d_ 

dz 



{z^ - If 



d 

dz 



d 



Hence, 



d 

dz 



z^ - lfF{z) 



2lil + l)iz'-l)Qiiz)—Qiiz 

dz 



d \ °° 

[z'-i)j-^\TX2i+i)mz)f 



Now recall that 



Y.{2i+i)mz)f 



1=0 



z^ - 1 



(Bl) 



(B2) 



(B3) 



(B4) 
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[This can be obtained by squaring both sides of Eq. ( |A18| ) and integrating over the variable 
t from —1 to 1.] By substituting Eq. (p^ ) in Eq. (|B3| ) , we obtain 



d_ 

dz 



[z^ - lfF{z) 
where we have used 



[z^-l 



I — 

dz 



Qo{z) 



1 4 V z-1 



1, z + l 

- In 

2 z-1 



2z , z + l 
h in 

2-1 z-1 



Integration of Eq. ( PSQ leads to 



z + l 



[z' - lyFiz) = z\n — y - 2 + C , 



(B5) 



(B6) 



(B7) 



where C is a constant. Recalling that Qi{z) z ' ^ as z ^ oo, we find that the left-hand 
side tends to zero in this limit. Therefore C = 0. Hence, 



F(z) 



{z^ - ly 



1 9 

z m 2 

z-1 



2Q,{z) 



(B8) 



APPENDIX C: ABSORPTION CROSS SECTION IN SPHERICAL POLAR 



COORDINATES: A DERIVATION OF EQ. (5.1) 



We consider the divergence-free vector eigenfunctions of the Laplacian on the {p + 1)- 
dimensional Euclidean space. That is, we examine the solutions of the following equations: 



and 



VBa = . 



(CI) 



(C2) 



Let X = (xi,X2, ■ ■ ■ iXp+i) be the Cartesian coordinates. There are plane wave solutions of 
the form 



(C3) 



with u > 0, where A; is a unit vector and k = uk, and e is a polarization unit vector 
orthogonal to k. We define 



for any vectors B^^^ and Bj^^ on this space. Then 

^^(C;^k)^^(C;^'fe'e'))^ ^ (27r)P+i| ■ e' 5^'+^ (k - k') 



(C4) 



(C5) 
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Now, suppose that B^''^'^^ are solutions of Eqs. (|Cl| ) and (|C^ ) with a discrete label k 
satisfying 

(^(P;..)^^(P;..'«'))^ = 27r5«,,5(^ - u') . (C6) 

Suppose that the plane wave solutions are expanded as 

^)^(P;..)^ (C7) 

re 

where a(a;, A;, e; k) are constants. As can be seen from the mode analysis in Section y, a 
solution of Eqs. ( piP and (|C2D satisfies the normalization condition ( p6|) if its incoming 



angular components at large r are Bi ~ r~^/^e~*'^''li, with the spherical harmonic on p- 
sphere Yi being normalized by / dVtpY^Yi = 1. (The component B^ can be shown to be 
suppressed by an extra power of r.) Thus, if we define the flux so that the solution B^''^^'^'> 
has a unit flux, then its content of the incoming mode B^^''^'^^ is \a{u!, k,e, k)\'^ per unit 
time. 

By using orthonormality of these solutions, we find 

a{uj,ke; K)27r6{uj-uj') = (^C^';-), ^(f^;-'^^'))^ . (C8) 
From this one readily finds 

Jo 

(C9) 

Now, the average of ejC^^'^'^eie'"''^''^'^' over the polarization vectors is 

-{5,1 - kfki)^^''^-^^-^'^ = -{5,1 - —M)e'-''^-^^-^'^ , 

p p UJ^ 

where the partial derivative d, is with respect to x whereas the partial derivative d[ is with 
respect to x'. Since the modes B\[' ^'^^ are divergence-free, these derivatives can be dropped 
when we average Eq. (|C9| ) over the polarization vectors. By averaging Eq. ( P9| ) over e and 
k and then using 

j du'J'Pdke'^''^-^''-'''^ = {27if5P-^\^ - x') , (CIO) 

we find \a{uj] the average of \a{uj, k, e; k)P over k and e, as 

\a{u;; k)\'u;^5{u - uj") = ^C^.-''^))^ . (cil) 

pilp 

Then from Eq. (|C^ ), we find 

Kc.;«:)p = ^. (C12) 
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Thus, if a plane wave has a unit flux, it contains the incoming mode with label k given by 
(C12) per unit time on average. 

The angular components Ai of physical modes I are gradients of scalar spherical har- 
monics. On the other hand, physical modes II are divergence-free vectors on the p-sphere. 
Let m/^"* and m/^^ be the multiplicities of these two types of modes with a given angular 
momentum /. Then the number of physical modes I with angular momentum / is m/^^ 
and that of physical modes II is . Then, if there is a plane wave with a unit flux, the 
content of physical modes I and II with angular momentum / are [{2ttY / pVLpUj'^]Ml^'^ and 
per unit time, respectively. Therefore, if the absorption probabilities of 
physical modes I and II by a spherical black hole are vl^"* and v'f'\ respectively, then the 
absorption cross section is given by Eq. (|5 
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